Avalanche dynamics driven by adaptive rewirings in complex networks 
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We introduce a toy model displaying the avalanche dynamics of failure in scale-free networks. In the model, 
the network growth is based on the Barabasi and Albert model and each node is assigned a capacity or tolerance, 
which is constant irrespective of node index. The degree of each node increases over time. When the degree 
of a node exceeds its capacity, it fails and each link connected to it is is rewired to other unconnected nodes 
by following the preferential attachment rule. Such a rewiring edge may trigger another failure. This dynamic 
process can occur successively, and it exhibits a self-organized critical behavior in which the avalanche size 
distribution follows a power law. The associated exponent is T ~ 2.6(1). The entire system breaks down when 
any rewired edges cannot locate target nodes: the time at which this occurs is referred to as the breaking time. 
We obtain the breaking time as a function of the capacity. Moreover, using extreme value statistics, we determine 
the distribution function of the breaking time. 

PACS numbers: 89.70.+C, 89.75.-k, 05.10.-a 



Complex systems are composed of many constituents that 
interact with each other or adapt to external perturbations 1 1 , 
0]. Recently, there have been increasing attempts to describe 
such systems in terms of networks |3, 4, 5, 6], where nodes 
and links represent constituents and their interactions, respec- 
tively. Many complex networks in real systems follow a 
power-law Pd{k) ^ k^^ in the degree distribution, where de- 
gree is the number of links connected to a given node |7]. 
Such a network is called scale-free (SF) networks. SF net- 
works are ubiquitous in real world, whose examples include 
the Internet and the world-wide web, the metabolic networks, 
the protein interaction networks, the co-authorship networks, 
etc. 

It was studied that complex networks are robust against the 
random removal of nodes, however, they are vulnerable to the 
intentional removal of nodes with high degree |0]. More se- 
vere damage can be caused by triggering a few nodes, but 
the failure propagates to other nodes in a cascading manner. 
Avalanche dynamics occurs frequently in complex networks 
due to the small-world feature of the system. The black- 
out of power-supply network in the United States in 1996 
and 2003 is a typical example of such a cascading failure 
in complex networks [9]. Internet traffic is another exam- 
ple. In October 1986, during the first documented Internet 
congestion collapse, the speed of the connection between the 
Lawrence Berkeley Laboratory and the University of Cali- 
fornia at Berkeley, which are separated by only 200 meters, 
decreased by a factor of 100 flu). These are only a few in- 
stances; many others can be found in various systems such as 
cultural fad, earthquake, etc. Since the subject of avalanche 
dynamics in complex systems is intere sting an d intricate, it 
has been studied extensively lUlllll flirM llSlliell. 

Here we introduce a toy model exhibiting the avalanche dy- 
namics in complex networks. In this model, when a node fails, 
the links connected to it are rewired to other nodes, thereby 
preserving the total number of links during the avalanche dy- 
namics. Due to the dynamic rule, the overload imposed on a 
specific node is shared among the other nodes globally as in 
the case of global load sharing in a fiber bundle model. 

Let us begin with the introduction of the model we consider 



here. The model is based on the Barabasi and Albert (BA) 
model 1 7], in which at each time step, a new node is added and 
its links are connected to m distinct existing nodes by follow- 
ing the so-called preferential attachment (PA) rule. The newly 
added node is connected to node / whose degree is fc, with a 
probability 11, = mki/Y.jkj. In this case, the number of nodes 
at f = is taken as m and they are fully connected to each 
other. Our model is modified from the BA model as follows: 
We assign a capacity or tolerance denoted by a to each node 
/. It is constant independent of the node index /. The capac- 
ity represents the maximum number of connections that can 
be sustained by a node. As time goes on, the degree of each 
node increases. When the degree of a node exceeds its capac- 
ity o, it is considered to be overloaded and is deleted from the 
system. Then all the links connected to the failed node are 
rewired to the other remaining nodes by applying the PA rule. 
By other remaining nodes, we mean the nodes in the system 
except the overloaded nodes. Multiple connections are not al- 
lowed. The rewiring process, shown in Fig. 1, may trigger the 
avalanche dynamics: when the degree of a node that receives a 
rewired link of an overloaded node exceeds a, i.e. it becomes 
O + 1, it fails and its links must be rewired again. This process 
repeats until all the overloaded nodes are eliminated. Thus 
the degree of each node that remains after the completion of 
the avalanche process does not exceed the capacity o. In the 
absence of a target node to which a rewiring link should be 
connected, the dynamics process is terminated and the entire 
system is considered to be collapsed, and the time at which 
occurs is referred to as the breaking time denoted by th. It is 
noteworthy that the avalanche process does not spread locally 
from the triggering node, but it occurs across the entire sys- 
tem. Moreover the number of links is preserved during this 
process since the load of the failed node is distributed to other 
nodes globally. Otherwise, the system breaks down. 

Using the model, we perform numerical simulations for 
various values of a. First we considered the breaking time as 
a function of the capacity o. As shown in Fig. 2, we observe 
that that the mean breaking time (f/,) averaged over different 
ensembles increases at a higher rate than the double exponen- 
tial function. This result implies that a small increment in the 



2 



9 




5 

FIG. 1: Illustration of the dynamic rule of the model for the case 
of a = 5. When node 9 is newly added and connected to node 1, 
the degree of node 1 exceeds its capacity a = 5. Then the node 1 is 
overloaded and deleted from the system, and the links i\ £(, (bold 
lines) connected to it are rewired. For example, the link £2 is rewired 
from 1 to one of the nodes 5, 6, 7, or 9. For example, it is rewired to 
node 6, pivoted on node 2. The target is selected according to the PA 
rule. 
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TABLE I: Summary table : capacity (a), the mean breaking time 
({th}), the root-mean-square of the breaking time (e) and the esti- 
mated values of the constants a, b, and r in the probability distribu- 
tion function F(f) = 1 — exp[— exp(a — fe^'')] 

capacity of each node significantly enhances the tolerance of 
the system. That is because the overload is shared globally. 
Since the breaking time increases very rapidly, it is difficult to 
perform numerical simulations for large values of a. 

The breaking time fluctuates considerably in Fig. 2, wherein 
the error bar increases with increasing a. Explicit values are 
listed in Table 1 . Thus, it would be meaningful to study the 
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FIG. 2: Plot of the logarithm of the mean breaking time (r^) as a func- 
tion of capacity a on semi-logarithmic scale. The linearly increasing 
behavior in this plot means ~ expexp(o). The data show that 
the mean breaking time seems to increase at a higher rate than the 
double exponential function. 

distribution of the breaking time since it is a relatively sensi- 
tive probe of the underlying breaking mechanism, which is of- 
ten used as a tool for reliability analysis. It is conventional to 
consider the probability that a network fails at a time f , or 
less. Then it is expressed as F{t) = fQPi,{th)dti,. In Fig. 3(a), 
we plot this distribution numerically, and observe that it fits 
well to the double exponential form, 

F(f) = 1 -exp[-exp(fl (1) 

which is one of the known functions in extreme value statis- 
tics. The constants a, b and r depend on o. Their val- 
ues are listed in Table I. The double-exponential functional 
form appears in mechanical failure problems in the duc- 
tile regime |17]. This functional form is different from the 
Weibull distribution that occurs in the brittle regime and fiber 
bundle model CI] 

Next, we study the avalanche size distribution in which the 
avalanche size represents the number of nodes that fail suc- 
cessively. This distribution follows the power-law, 

Pa{s)-s-\ (2) 

with the exponent T « 2.6(1) (Fig.|^. The exponent value is 
close to that obtained from the fiber bundle model. We deter- 
mine the mean avalanche size as a function of triggering time. 
As shown in Fig. 5, the mean avalanche size is independent of 
the triggering time except in the small t regime. 

The exponent of the avalanche size distribution of the toy 
model is close to the one obtained in the fiber bundle model 
in SF networks |19]. In the original one-dimensional fiber- 
bundle model, an external force is evenly distributed to all 
the nodes in the network. For O; of node /, a threshold value 
is assigned against failure. When the load is larger than the 
threshold, the node fails irreversibly and its load is equally 
distributed to its remaining nearest neighbors. This process 
continues until no failure occurs. Thereafter the process is 
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FIG. 3: (color online) (a) Breaking time distribution Pb(th) for CJ=13 
on linear scales, (b) Probability of failure at time / or less, and F(t) 
as a function of /. Inset: Determining whether F(t) fits to the double 
exponential function. 
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FIG. 4: Avalanche size distribution /'^(i) for various o = 12 (□), 
13 (O), 14 (A), 15 (V), and 16 (O) on double logarithmic scales. 
The slope of the linear fit (solid line) is —2.6. The data are averaged 
over 10^ configurations for a = 12, 13, 14 and 15, but over 55000 
configurations for o = 16. 



FIG. 5: Average avalanche size {s) as a function of time for ct = 16. 



repeated by applying a bigger force to the remaining nodes. 
The study of the fiber bundle model has been extended to 
the case of complex networks. When a critical force is ap- 
plied, the avalanche size distribution follows a power law. Re- 
cently Kim etal. 1 1 8 ] investigated the fiber bundle problem by 
considering various complex networks such as the Erdos and 
Renyi random network 1 20 ] , the small-world network 1 2 1 1 and 
the SF network. Their result indicates that the patterns of the 
avalanche dynamics occurring in such complex networks are 
almost the same due to the small- world property: The effect of 
the local load sharing rule is negligible in complex networks. 
Therefore the exponent of the avalanche size distribution re- 
duces to the mean field value, i.e., T = 5/2. 

As shown in Fig. 5, the avalanche size is independent of 
time. This suggests ineffective correlation between the de- 
grees and the entire system is considered to be homogeneous. 
Moreover, the avalanche dynamics proceeds globally by the 
rewiring dynamics. Thus, the exponent of the avalanche size 
distribution of the toy model reduces to the mean field value. 
Although the toy model has not been applied to real world 
systems, the dynamic rule of the rewiring edges in the model 
is rather unique; it may reflect the adaptive behavior of each 
individual (node) in the event of a failure. In such a case, the 
adaptive behavior does not relieve itself, rather it may lead 
to other successive failures. This phenomenon can often be 
observed in complex systems. For example, when a city en- 
counters a shortage of electric power, the cuiTent power grid 
system is designed to compensate instantly by drawing power 
from neighboring cities. This may cause another blackout or 
a cascading blackout throughout the country. In order to pre- 
vent such avalanche dynamics induced by the adaptive activ- 
ity, failure should be localized. 
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